Covariant kinematics and gravitational bounce in Finsler space-times 
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The similarity between Finsler and Riemann geometry is an intriguing starting point to extend 
general relativity. The lack of quadratic restriction over the line element (color) naturally generalize 
the Riemannian case and breaks the local and global symmetries of general relativity. We investigate 
the covariant kinematics of a medium formed by a time-like congruence. After a brief view in 
the general case we impose particular geometric restrictions to get some analytic insight. Central 
role to our analysis plays the Lie derivative where even in case of irrotational Killing vectors the 
bundle still deforms. We demonstrate an example of an isotropic and exponentially expanding cross- 
section that finally deflates or forms a caustic. Furthermore, using the 1+3 covariant formalism we 
investigate the expansion dynamics of the congruence. For certain geometric restrictions we retrieve 
the Raychaudhuri equation where a color-curvature coupling is revealed. The condition to prevent 
the focusing of neighboring particles is given and is more likely to fulfilled in highly curved regions. 
Then, we introduce the Levi-Civita connection for the osculating Riemannian metric and develop a 
(spatially) isotropic and homogeneous dust-like model with a non-singular bounce. 

Keywords: Finsler, Raychaudhuri, gravitational collapse, cosmology, non-singular, bounce 



I. INTRODUCTION 

In Einstein's theory of gravity the 4-dimensional space- 
time is described by a pseudo-Riemannian manifold. Lo- 
cally, the curvature effect is diminished and the Lorentz 
symmetry is always restored. However, when we investi- 
gate long range phenomena, the distance module is non- 
invariant under the 5*0(1, 3) group since we depart from 
the Minkowski line-element. The most relevant manifes- 
tation of this Lorentz symmetry breaking is the gravita- 
tional redshift of light. When we consider hypothetical 
geometric modifications of general relativity (GR) we can 
violate the Lorentz symmetry locally and/or modify the 
way we break it globally. These scenarios are the com- 
mon playground for phenomenological models of quan- 
tum/emergent gravity. 

From a geometric perspective, we can introduce such 
GR modifications by dropping the restriction that the 
line-element should only depend on quadratic terms with 
respect to the coordinate increments. This non-trivial 
"symmetry breaking" is the essence of the velocity de- 
pendent geometry named as Finsler [H-Q. In particular, 
the lack of the quadratic restriction directly introduces 
Lorentz violations (LV) [IHl] since the space-time locally 
is no longer invariant under the boosts of special rela- 
tivity. Also, the global characteristics of the space-time 
geometry are modified as a consequence of the deformed 
length interval. 

Extending the notion of distance in the field of non- 
quadratic metric functions brings into play an extra ge- 
ometric property, apart from the curvature, related to 



departures from the quadratic measurement. This extra 
property of the Finsler manifold, often referred as color 
0, Q > directly introduces LV effects within a geometric 
framework. It may be considered as an effective mani- 
festation of quantum/emergent gravitational phenomena. 
Moreover, the color of the manifold enforces us to intro- 
duce two types of covariant derivatives in order to prop- 
erly define the parallel displacement. The first covariant 
derivative (derivative of the horizontal space) generalizes 
the Riemannian derivative, while the second (derivative 
of the vertical space) is a pure Finslerian entity. Ap- 
parently, this "exotic" differentiation of Finsler geome- 
try gives rise to three distinct Ricci identities that result 
from mixing up the vertical and horizontal derivatives. 
This set of covariant identities encodes all the informa- 
tion about the kinematics of Finslerian flows. 

Finsler spaces can be equivalently described as the ge- 
ometry of a velocity dependent metric where the distance 
between two points is invariant under reparametrization 
of the connecting curve. In the literature velocity depen- 
dent metrics have been studied from the onset of gen- 
eral relativity, especially in geometric optics and ana- 
logue gravity. One of the main "privileges" of Finsler 
space-times is that they naturally induce modified dis- 
persion relations [ll| . Recently, this type of geometry ap- 
peared in different perspectives of quantum gravity (QG) 
(see for example [10| and references therein). For exam- 
ple: in Horava-Lifshitz gravity rays may follow Finsler 
geodesies [12fl, within the context of the the "stringy" 
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space-time foam where the D-particles recoil with the 
world-sheet in Cohen and Glashow's very special rel- 
ativity , in bi-metric theories of gravity [l|J EH and in 
holographic fluids [l6[ . There are also studies concerning 
extensions of relativity theory within the framework of 
pseudo-Finsler structures dating back to the 40 's [T7l.[l8j 
with several issues still remain open [l9j . 
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In GR the internal motion of a flow is "encrypted" in 
the failure of the deviation vector to be parallel trans- 
ported along the congruence. Following the mathemati- 
cal definition of the Lie drag along the flow we can track 
the evolution of expansion, shear and vorticity with the 
aid of the Ricci identity. This elegant and straightforward 
derivation leads to Raychaudhuri's equation that played 
a keynote role in the development of GR (for some re- 
cent works see [2l|). It shed light to crucial questions 



like gravitational collapse and singularities, gravitational 
redshift and lensing, accelerating expansion and inflation. 
This discussion is carried on to various branches of QG 
where intriguing effects are reported [22j]. In the same 
line, the recent entry of Finsler geometry in emergent as- 
pects of QG creates the need for further insight in the 
volume evolution of a congruence. 

Motivated by the mentioned research works we will 
investigate the deformation of time-like congruences fol- 
lowing E. Cartan's metrical approach [l|, |5|, l2C§ . The kine- 
matics of Finsler geometry includes several complications 
coming from the three different Ricci identities and the 
anisotropic corrections in the Lie and absolute derivative. 
One must solve the global problem of neighboring curves 
and their deviation vector by imposing specific symme- 
tries for the medium. In particular restrictive cases we 
demonstrate how we can retrieve a closed form expression 
for the expansion of neighboring particles. The general- 
ized Raychaudhuri's equation is derived, where a color- 
curvature coupling is revealed that becomes more im- 
portant in highly curved regions. The simple case of a 
dust-like fluid moving along irrotational geodesies gives 
back a "counter gravity" condition where the focusing of 
matter is prevented. Finally, the limit of osculating Rie- 
mannian metric is considered for a Randers space-time. 
In highly spatially curved regions the Finsler contribu- 
tion can dominate and lead to a bouncing Friedmann- 
Robcrtson- Walker (FRW) model. 



II. FINSLER SPACE-TIMES 

Let <t(t) be a smooth curve in a manifold M. connect- 
ing two space-time points. The distance between between 
the points is given by the following integral over the con- 
necting curve and its tangent direction 



£= / F(a(T),a(T))dr. 



(1) 



The metric function F is defined over the tangent bundle 
TM- By imposing that the distance ([!} is independent 
of the curves parametrization r, we get 



F(x,Xy)=XF(x,y), A > 0. 



(2) 



where y a = is the tangent vector to the cr(r)-curve. 
In other words, the metric function F(x,y) is homoge- 
neous of degree one with respect to the coordinate incre- 



ments y a . Using Euler's theorem of homogenous func- 
tions we retrieve the formula 



F 2 (x,y) = g ab (x,y)y a y b , 



(3) 



where g ab is the metric tensor and is given by the relation 



9ab{x,y) 



1 d 2 F 2 



2 dy a dy b ' 



(4) 



Furthermore, a Riemann manifold is retrieved if F(x,y) 
is quadratic with respect to y a . Then, the double par- 
tial derivative of the previous relation removes the in- 
dependence of the metric. In that sense, Finsler ge- 
ometry generalize the Riemannian case within the field 
of non-quadratic metric functions. The Euler-Lagrange 
equations for the extremal of F(x, y) give back the usual 
relation 



dr 



2G a = 0, 



(5) 



where the second term in the lhs is usually called the 
spray induced by F(x, y) and is related to the Christoffel 
symbols by the relation 



(6) 



The extremal curves © define a space of paths for the 
Finsler manifold. 

The essential difference in the Finslerian approach is 
that the quadratic restriction on the metric function is 
dropped. Consequently, the metric tensor depends on 
the variable y a and each tangent space T X M is no- longer 
equipped with ellipsoidal unit balls. Instead, there will be 
a locus given by the restriction F(x,y) — 1 (see rel.([3])). 
This non-ellipticity of the unit balls gives an extra prop- 
erty to the Finslerian manifold, apart from curvature, 
referred as color. The main quantity that measures the 
color of our structure is the Cartan torsion tensor 



Cab 



dg a b 
dy c 



(7) 



since the condition C a bc = implies a Riemann space. In 
that sense, a Finsler space is a colorful curved manifold, 
while a Riemann space is curved but entirely white [361 ]. 
Recently, this extra property of the J- n manifold drew 
some attention in the community of theoretical physics 
since it is related to Lorentz violating scenarios (see for 
example 0], || and fllj|). 

We consider the position space together with the y a 
increments, that solve relation ([5]), as the fundamental 
variables of the physical space-time. In other words we 
replace the Riemann space-time with a Finsler one where 
the position space gives place to the element of support 
( x ,y) p|.[20|. Thus, tensor fields will depend on the co- 
ordinates of the tangent bundle TM. induced by the local 
chart of the base manifold A4. In the spirit of [Hj|, the 
square magnitude of a first rank tensor u a in x is given by 
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g a b{x,y)u a u b and it separates to three classes, spacelike, 
null and timelike for negative, zero and positive values 
respectively. 

The tangent bundle TM. has a local coordinate base 
{gfs, 5^3"}- However, the elements do not trans- 
form properly under coordinate transformations on M.. 
This problem is solved by using instead the adapted frame 
{^>^}»with 



6 

Sx a 



d 
dx a 



dy a '' 



(8) 



where N„ 



9G° 
dy" 



is the non-linear connection. Under 



the same argument, the adapted frame in the cotangent 
bundle T* M is {dx a : 5y a }, with 



6y a = dy a + N a b dx a . 



(9) 



The presence of the non-linear connection imposes that 
the adapted frame is non-holonomic in the sense that the 
following commutation formulas hold 



[— -] 

l &r a ' Sx bi 

[— —) 

L fe a ' dy bl 

d d 



R c 
G 
= 0, 



d 

ab dy~ c ' 
d 

ab dx c 



where we used the definitions 

5N\ 5N" 



pa 

n be ~ 



Sx c 



Sx b ' 



and 



G" 



be 



dy c 



(10) 

(11) 
(12) 

(13) 
(14) 



Note, that the tensor field (TT3"| is also present in the 
Riemann limit and reflects the gravitational field. On 
the other hand, the G a b components monitor the "mess 
up" between the coordinates of the base manifold A4 and 
the elements of the tangent bundle y a . 



III. ASPECTS OF FINSLERIAN 
CONGRUENCES 

Consider a family of curves "f(s) and their tangent vec- 
tor field u a in a manifold M.. The parallel displacement 
of a vector X a along the w a -bundle is given by the fol- 
lowing formula 



DX a dX a 



ds 



ds 



pa 



(15) 



The components P a are to be specified by the geometric 
properties of M.. We consider the mathematical defini- 
tion of the deviation vector £ a as the tangent of a second 



flow given by the Lie drag of a connecting curve along the 
congruence. The failure to parallel transport the devia- 
tion vector along the flow monitors the internal motion 
of the u a -family, namely 



Dt] a 
ds 



(16) 



The tensor fields in the rhs of the above can take explicit 
expressions depending on the geometric assumptions we 
impose. The common example is the Riemann case given 
in the next subsection. 



A. White curved space-times 

When the manifold is Riemannian the parallel dis- 
placement (1151) involves only corrections given by the lin- 
ear connection. Taking into account that the deviation 
vector is Lie transported along the flow, £ u £, a = 0, the 
absolute derivative of £ a can be treated as 



ds 



e = £, b v b u a . 



(17) 



Combining the previous expression with relation (|16j) we 
find that only first order terms are involved in our ap- 
proximation. Using the 1+3 covariant formalism (37| . we 
can write into its irreducible parts the second rank tensor 
that encodes all the relevant information 

Bab = Vfciia = ^<dh ab + a ah + co ab + A a u b , (18) 

where we use the definitions with respect to the 3-D 
covariant operator D a = h a b W b , for the expansion of 
a volume element = D a u a , for the shear a ab = 
2T)( b u a ) — ^Qh ab , for the vorticity uj ab — D[ h it a ] and the 
four acceleration A a = u a [23| . 

In Riemann geometry the covariant derivative on a vec- 
tor field fails to commute due to the curvature of space- 
time. The tensorial Ricci identity reads 



2V[ Q V b ]U c = RabcdU 



(19) 



and decomposes to three propagation and three constrain 
equations [23j for the irreducible parts of relation (|18|) . 
The most relevant equation of this set is Raychaudhuri's 
equation that monitors the evolution of the average vol- 
ume of a spatial 3-D element dV. It results directly from 
the Ricci identity flT9"|) by taking the trace and projecting 
along the w a -family 

<d + \o 2 = -R ab u a u b -2{<r 2 -uj 2 ) + D a A a -A a A a , (20) 

where a 2 — a ab a ab /2 and w 2 = w ab uj ab /2 and the last 
two terms are representatives of external forces. Given an 
initially expanding phase of the flow (9 > 0) we directly 
conclude that positive terms in the rhs will accelerate the 
expansion and negative ones will decelerate. The state- 
ment is reversed in case of a contracting phase. Notice, 



4 



that the evolution law (|20j) is a pure geometric relation. 
It is valid for any theory of gravity that assumes a Ric- 
mann manifold for the position space. We can indirectly 
connect the evolution formula (f2"U|) with the energy and 
momentum of matter through the field equations and the 
conservation relations. In the context of GR Raychaud- 
huri's equation reflects the attracting nature of gravity, 
since for ordinary matter the curvature contribution al- 
ways assists the contraction [21j . 



B. Parallel displacement and Ricci identities 

As we discussed, in Finsler spaces tensor fields no 
longer depend only on the position but instead they de- 
pend on the element of support (x,y) [H, Therefore, 
the properties of the tangent bundle TM is of central 
importance. The absolute differentiation of a vector field 
X a will contain extra contribution along the tangent dis- 
placement dy a . In a local expression we can write down 



DX a = dX a + T a bc X b dx c + C a bc X b dy c 



(21) 



where the coefficients C a bc and r a 6c are functions of the 
element of support. We will restrict our analysis only 
along the normalized time-like direction u a = This 
congruence plays an analogous role with the fundamental 
observer of the 1+3 covariant formalism. 

By imposing the metricity condition Dg ab = 0, along 
any direction u a we can recast relation (|2T|) in a covariant 
form by using the horizontal and vertical split 



<ra _ Dxa va 
X = ~~ds~ — I' 



where 



X a , 
X a \ b 



sx a 

dX a 



dy b 



X a \ b F 



s~ia vc 



Dl b 
ds 



(22) 

(23) 
(24) 



L a bc l b dx c . 
implies l a ^ b = and 



and l a = y a /F is the unit vector in the direction of the 
element of support (x, y) while Dl a — dl a 
Note, that for X a = nl a relation i 

the spray coefficients ([5]) determine the integral curves. 
We will refer to the operator (1231) as the horizontal co- 
variant derivative and to (|24p as the vertical one. Also, 
the connection coefficients in (|24l) satisfy the following 
relation 



ja _ 1 a 
L bc - 77.9 



Sgdb . Sg dc 



8x c 



5x b 



5gbc 
5x d 



-Er ■ (25) 



In that case a parallel displaced vector field DX a = 
keeps its length invariant. The covariant expression (|2"2")l 
is the evolution law for any first rank tensor. The vertical 
part (|2~4"|) is a direct result of the distorted ellipsoid con- 
dition ([3]). Together with the non-linear connection © 
monitors the LV effect on the evolution of the X"-flow. 



The splitting of the propagation relation (|22j) allows 
us to formulate a curvature theory in the footsteps of 
Riemann geometry. Using the vertical and horizontal 
covariant derivatives we retrieve the following generalized 
Ricci identities 



va va vdjj a va\ r>d 

A iu„ A i j. — A K d bc — A \dH 



\b\c 



x a \ 



\c\b 



bc > 



vd q a 
A D d be! 



(26) 
(27) 



and 



A | b | c -A c | b - A r dbc -Ji | d O bc -A \d-r bc , (*o) 

where the torsion-like tensor field P a bc is given by the 
relation 



pa _dN% 



dy c 



L a 



cb- 



(29) 



The Ricci identities (|26p - ([2T|) enclose the relevant infor- 
mation about the non-linear kinematics of the integral 
curves since the torsion-like tensors C a bc , R a bc and P a bc 
are fully determined by the fundamental function F(x, y). 

The irreducible decomposition of relations (|26p - ([28|) 
will give back a set of covariant relations [30(. However, 
this is not enough to determine the actual deformation of 
the X a -curves. One must further proceed to connect the 
aforementioned relations with the actual internal motion 
of the congruence. In the general case, it seems difficult 
to retrieve closed form expressions for the evolution of 
expansion, shear and vorticity. Therefore, to get the de- 
formation of a cross-section we need to solve the global 
problem and track back the kinematics of the medium, 
either by considerable simplifications or by using numer- 
ical techniques. 

To further dig into this problem, consider a family of 
time-like curves and their tangent vector u a in a Finsler 
space-time. By virtue of the absolute differentiation (|2"2")l 
and and using (|16[) , the propagation equation for the de- 
viation vector is written in the form 



nib 

c lb u b +e\ b F— = B%e 



A abc ee 



(30) 



In order to get further analytic insight we must correlate 
the distortion tensors B ab and A abc with the 4-velocity 
of the ^-congruence. Depending on the definition of 
the deviation vector we can drop or kee p h igher order 
terms in the rhs of the previous relation [24j |. However, 
in the lhs the vertical derivative involves such terms in the 
kinematics. This fact leads us to the usual mathematical 
definition of the deviation vector as the Lie transported 
field along the relevant congruence. 



IV. LIE DERIVATIVE AND KILLING VECTORS 

The most closely related to the Riemannian case defi- 
nition of the Lie derivative results from the infinitesimal 
transformation 



x a =x a + v a {x)dT, 



(31) 
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where dr is an infinitesimal constant The above dis- 
placement assigns at each point x a a shift on the base 
manifold dx a = v a (x)dr. Following the usual procedure 
for a vector field X a (x, y) we lead to the formula 

f) X a 

£ v X a = X a lb v» - v a ]b X b + ^r(^V). (32) 

The third term on the rhs of the previous relation comes 
from the preferred direction imposed by the element of 
support. However, the reader should keep in mind that 
relation (|32"T) is a special case, since the Lie derivative 
is defined along the vector v a (x) that is independent of 
the coordinate increments y a . The condition £ v g a b = 
implies the equation for the Killing vectors 

v alb + v bla + 2FC abc v c {d l d = 0. (33) 

Thus, the Killing solutions are kinematically enriched 
compared to general relativity since Cartan's torsion ten- 
sor Cabc and the element of support implies expansion, 
shear and vorticity for the bundle. 

Consider the 4-velocity of a time-like congruence u a 
together with the deviation vector £ a and let us assume 
that they both depend only on the position x a . In that 
case, relation (|32l) takes the simplified form 

e lb u b - u a lb e = o. (34) 

Then, for a parallel displaced element of support, — 
0, we can drop the second order terms in relation (1301) . 
Assuming a geodesic time-like congruence and using (|30j) 
together with ([33]) and (|34f we get 

e a6 + FC abc {a c d + uj c d ) i d = 0, (35) 

where we define C abc — h a d h b e hj Cdef and & ab — 
■|0/i a b + <7 a b in analogy to the Riemann case for the 
horizontal covariant derivative, while £ a — h a b l b is the 
space- like part of the unit vector l a . Thus, for a time- like 
congruence u a (x) that is a Killing vector and the devia- 
tion vector depends solely in the position , the expansion 
and shear are non-zero if the element of support points 
to a preferred direction in space. Note, that shear free 
and irrotational Killing vectors of the previous example 
are always non-expanding. 

V. FINSLERIAN CURVES AND VOLUME 
EVOLUTION 

The only geometric tool in hand to measure the rel- 
ative motion of a family of curves j(s) is the parallel 
displacement of the deviation vector. As we already 
mentioned this displacement is directly related to the co- 
variant derivative of the tangent vector u a . However, in 
Finsler geometry the covariant derivative is replaced by 
the horizontal and vertical decomposition of the absolute 
derivative. This is a natural result, since the abandon- 
ment of the quadratic restriction forces us to replace the 



position space x a with the element of support (x a ,y a ). 
Therefore, to isolate the volume evolution of a Finslerian 
congruence we need to take into account all the curva- 
ture together with the torsion tensors (j2"6"|) - (f25|) and up 
to second order terms in relation (fT6|) . The analytical 
treatment of this problem seems a very complicated task 
but we can get some insight through some geometric as- 
sumptions. 

Firstly, let us prove that second order terms in relation 
(fT6|) enter the kinematics. Consider a time-like flow u a 
and its deviation vector £ a that depends only on x a . The 
parallel displacement of the unit direction l a along u a 
decomposes to a component along and perpendicular to 
the deviation vector £ a 

Dl a 

— =r + f a . (36) 

To focus on the second order terms we assume that the 
transverse part of (|36|) vanishes, f a = 0. Then, by virtue 
of (US]),® and dM|) we get the following relation 

^ = u a lb e + fFC a bc e^ c , (37) 

which proves our previous statement that second order 
terms are naturally involved in the evolution of the de- 
viation vector. It is evident from relation (|37p that an 
analytic expression for the expansion, shear and vorticity 
of the congruence is unapproachable. The second order 
terms of £ a complicate the picture and even in simple 
geometric configurations non-linear effects will emerge. 
Even if we consider the isotropic condition for the first 
order terms u a \ b = Bh ab , Cartan's tensor will change 
the initial direction of £ a . Therefore, an isotropic cross- 
section will eventually be distorted. 

A. Spherical symmetric evolution 

Assume that Cartan's torsion tensor takes its reducible 
form H 

Cabc = C(at)bc), (38) 

where C a = C a bc g bc is often referred as Cartan's torsion 
vector. Also, \) a b = 9ab + cl a h is the angular tensor that 
projects orthogonal to y a and c = =pl for a space- like 
or time-like unit direction respectively. Apart from the 
anisotropic contribution from Cartan's tensor in relation 
(|3T|) there is an isotropic effect along the direction of £ a . 
In particular, the relevant isotropic evolution is when C a 
is parallel to the deviation vector C a — C£ a and the 
simple relation u a \ b = Bh ab holds. Then, using relation 
(|38p we arrive at the expression 

-^- = (B + Ce)il a , (39) 

with C — fFC. Thus, in the above example the deviation 
vector always points at the same direction as it is parallel 
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FIG. 1: Free scale volume evolution for an initially expanding 
cross-section with respect to the parameter of the fundamen- 
tal tt a -congruence. The solid line corresponds to negative val- 
ues of C. An initially exponential expansion gradually turns 
to a decelerating phase that finally "freeze" to a constant 
value Vm = ( — §) 3//2 - On the other hand, positive values of 
C (dashed line) set an upper value s, for the arc-length where 
the volume element reach the caustic. 

displaced with respect to the 7(3) curves. An initially 
spherical symmetric region of space will sustain its shape 
as it moves along the it a -congrucncc. In that case, the 
magnitude of the deviation vector is always the radius 
r of an isotropic 3D-section. Contracting relation 
along £ a we derive the evolution law for the radius 

r = {B + Cr 2 )r, (40) 

where the dot operator denotes the absolute derivative 
with respect to s. The second term in the rhs of relation 
(|40jl involves non-linear corrections to the expansion dy- 
namics through Cartan's torsion vector. When Cartan's 
tensor vanishes, we retrieve the vorticity and shear free 
expansion of the Riemannian limit where B — ^0. In 
that case a constant expansion rate guarantees an expo- 
nential solution for the volume element. The exponential 
behavior is closely related to inflationary scenarios of the 
early universe. 

To demonstrate how a simple inflationary model can 
be affected by the non-linear Finsler contribution let 
B = const and C = const. This limit retrieves the afore- 
mentioned behavior for C — > 0. The solution of relation 
(|40p leads to the formula for the volume of an isotropic 
cross-section of the u°-congruence 

/ B \3/2 

y= U«--c ) ■ < 41 > 

where c± is a constant of integration. In the Riemannian 
approximation, the case of B < corresponds to an ex- 
ponential contraction which tends to the singular value, 
V — > 0. For the Finslerian volume element (|4"Tj) given 
the conditions B < and C > we get an additional 
phase of decelerating contraction for small values of the 
parameter s. On the other hand, for B < and C < 
the element monotonically contracts. In both cases of 



contraction the singular fate, V — >■ 0, is always unavoid- 
able. The inflationary behavior in general relativity is 
retrieved for B > where a 3-sphere eternally expands 
exponentially In relation (|4"Tj) , for B > and C < 
an initially inflationary expansion turns to a decelerat- 
ing one that leads to a "frozen" state, V m — (— ^) 3 ^ 2 - 
Furthermore, for B > and C > and restricting the 
values of V to be real, the parameter s is bounded by an 
upper value s* = -g In ^j- . The volume element tends to 
a caustic as we approach s* since its radius reach infinity 
within a finite time interval. The two cases of expansion 
are depicted in figurc-[TJ 

B. Raychaudhuri's equation 

The intrinsic non-linear nature of the Finsler manifold 
disapproves the extraction of the propagation equations 
for the kinematical properties of the u a -bundle. In the 
general case the evolution of the expansion cannot be 
expressed in closed form. This is evident even in the 
simplified example of the previous section (see relation 
pUf). The main suspect for this complication is the ver- 
tical part of the absolute derivative. In order to derive 
Raychaudhuri's equation we are forced to impose some 
geometric constraints. The first drastic assumption taken 
into account concerns the absolute derivative of the unit 
direction l a . 

Precisely, let us consider a family of time- like flow lines 
along which l a is parallel displaced, ^i— = 0. In that case 
the horizontal derivative along the congruence gives back 
the evolution of the deviation vector, £ a = £, a \bU b . The 
unit clement is written into its irreducible parts as 

l a =lu a +£ a , (42) 

and since the absolute derivative along u a vanishes we 
retrieve the expression for the 4-acceleration of the con- 
gruence 

lA a = -i {a} . (43) 

Note, that if the unit direction is purely space-like or 
time-like then the u a -flow is a geodesic congruence. Us- 
ing relation (|4"5|) together with the limit that the connect- 
ing vector £ a depends only on x, the Lie derivative (|54"]) 
implies that only first order terms are involved in the rhs 
of (|30|) . Moreover, we can track the internal motion of 
the medium with the aid of the horizontal derivative 

B ab =u a \ b . (44) 

The irreducible parts of the above expression monitor the 
expansion, shear and vorticity of the congruence. Their 
definitions are similar to the Riemannian case (fT5|) al- 
though the Levi-Civita connection V a is replaced by the 
horizontal operator. 

Furthermore, the time-like part of the contracted Ricci 
identity gives back the propagation equation for the ex- 



7 



pansion 



-RabU a U 

1 

1 



2(a 2 

-'<»>) 



u 2 ) - 1 ab u a u b 



(45) 



where a 2 = a ab a ab /2 and uJ 2 = uj ah uj ab /2 are respectively 
the scalar square magnitudes of the shear and vortic- 
ity, while 1 ab = C ac d l e R e c db represents a color-curvature 
coupling coming from the anholonomy of the basis and 
the intrinsic anisotropy of the metric. The last three 
terms are due to the acceleration given in relation (|43[) . 
In the rhs of relation (|45|) positive terms resist the con- 
traction of the time-like bundle while negative terms as- 
sist the focusing of the congruence. On the other hand, 
in an expanding phase positive terms accelerate the devi- 
ation of neighboring flow-lines while negative terms con- 
tribute with a decelerating effect. 



VI. COLLAPSE OF A DUST-LIKE FLUID 

In general relativity Raychaudhuri's equation is a nec- 
essary tool to study the behavior of a self-gravitating 
fluid. It has a purely geometrical origin although it is in- 
directly correlated to the energy and momentum of mat- 
ter through Einstein's field equations. When we consider 
irrotational (Lu a b = 0) perfect fluid configurations of ordi- 
nary matter the congruence always forms a caustic. On 
the other hand, if we replace the Riemann manifold with 
a Finslerian we are unable to derive in the general case a 
closed form expression for the expansion. Also, the rela- 
tion between matter and geometry remains in the sphere 
of speculation but with remarkable progress been made 

The derivation of relation (|45|) requires a parallel dis- 



placed element of support t~ = along the unbundle. 
Under this assumption the relevant Finsler entities that 
enter Raychaudhuri's formula are the Riemann curvature 
(l26f together with the two torsions (C a bc , R a bc ) and the 
4-acceleration. Furthermore, non-geodesic motion is nec- 
essarily related with the unit preferred direction through 
(|4"3"| . Since, for the parallel displacement (|2"21 the only 
operator is the horizontal differentiation (|23|) it is identi- 
fied as the covariant derivative. If we impose energy and 
momentum conservation for our setup, then we get 



rpab 



0. 



(46) 



with T ab representing an imperfect fluid. One can search 
in the identities of Finsler geometry for a second rank ten- 
sor that satisfies relation (l46l) . Keeping close to GR, the 
Bianchi identities for the Riemann curvature (|26|) with 
respect to our covariant differentiation are 



1?„ 



cd\e 



p b 

•Fa cf 



R„ 



R'.P. 



de\c 



cd a ef 



RL 



p b — n 



(47) 



The nearest curvature theory to general relativity is 
achieved if the space-time has the property P a bcd = 0. 
Then, after some short calculations we arrive at the fol- 
lowing conservation law 



R (ab) + R [ab] _ l R al 

2 3 



0, 



(48) 



where R ab — g Racbd and R = g R ab is the Ricci tensor 
and scalar respectively. Also, due to the torsion contri- 
bution in the curvature ([26]) the Ricci tensor decomposes 
into a symmetric and antisymmetric part. 

The most trivial choice is to assume that the symmetric 
and antisymmetric parts of the twice contracted Bianchi 
identities (|48l) are explicitly conserved. Consequently, 
relations (|46|) and (|48| lead to the field equations 



R 



(ab) 



-^Rgab 



i^T ab , 



(49) 



where k is the gravitational coupling constant, together 
with the constrain relation 



rt | fe — u. 



(50) 



The above field equations represent a subset of the gen- 
eral Finsler problem [l9|,[25j]. Although, they are compat- 
ible with the geometric assumptions we made to derive 
the evolution equation for the expansion (|45|) . Relative to 
the fundamental observer the energy-momentum tensor 
of a perfect dust-like fluid is 



T, 



ab = PU a U b , 



(51) 



with p = T ab u a u h standing for the energy density of mat- 
ter. Notice, that relation (|49p implies that the energy 
density will depend on the element of support. 

The conservation law (|46|) forces the matter distribu- 
tion (|5T|) to move along geodesies. Therefore, the 4- 
acceleration vector (|43l) should vanish. This implies that 
either the unit direction l a is purely space-like or time- 
like, or both its space-like and time-like parts are con- 
stant along u a . Then, by virtue of relations (|49| . (fSTj) and 
assuming that matter falls along irrotational geodesies, 
we recast Raychaudhuri's equation (|451) in the following 
form 



e + \e 2 



i 

'2 Kp - 



2a 2 - %abu a u b 



(52) 



The sign in the rhs of (l52"j) determines the condition for 
the formation of a caustic. In particular, the focusing 
of the world-lines is prevented only if the rhs of relation 
(|52p is dominated by positive terms, namely 



1 

2 KP 



2a 2 < -% ab u a u b 



(53) 



The term that can resist the gravitational pull of a 
geodesic congruence comes from the coupling between 
the Cartan torsion tensor and curvature. It fulfils the re- 
quirements only if % ab u a u b < 0. On the other hand, the 
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color-curvature effect adds to the self attraction of the 
medium when 1 a i,u a u b > 0. The more we deviate from 
the quadratic restriction the more important becomes 
the Cartan torsion and the coupling term gets stronger. 
However, even for slightly colored regions (C a b c <C) the 
prevention of the caustic becomes plausible in highly 
curved conditions since the Riemann tensor enters the 
rhs of relation 



VII. OSCULATING RIEMANNIAN FRW 
GEOMETRY 

Another method to explore some aspects of Finsler ge- 
ometry is the osculating Riemannian approach. This lim- 
iting process restricts our analysis to a region of F n where 
the coordinate increments y a depend on the position x a 
(see for example [l], I2H ] and for the kinematics of de- 
formable media j29[). In that case relation Q defines a 
purely Riemannian metric, namely 



r a b{x) = g a b{x,y(x)), 



(54) 



and the imposed metricity implies the Levi-Civita paral- 
lelism. Then, in our subregion we can use the Einstein 
field equations together with the energy-momentum con- 
servation by following the usual arguments 

G ab {x,y{x)) = K T ab (x,y(x)) , V b T ab {x, y{x)) = 0, 

(55) 

where G a b — Rab — \Rr a b stands for the Einstein tensor, 
while the nabla operator is the linear connection induced 
by the metric 

Let us consider the interesting class of Randers spaces 



F = a 



(56) 



with a = y a a by a y b standing for a pseudo-Riemannian 
metric, while (3 = b a y a is the Finslerian contribu- 
tion. The vector field b a introduces a preferred direction 
in space-time as a phenomenological consequence of a 
Lorentz invariance violation. The metric function (|56[) in- 
terfaces a Riemannian space-time with Finsler geometry 
in a simple manner. The study of such space-times gives 
many new possibilities since it monitors departures from 
specific Riemannian examples. Concerning the process 
of gravitational collapse in general relativity the most 
characteristic example is given by a spherical symmetric 
and isotropic dust fluid. In that case, the geometry of 
space-time is well described by the FRW metric 



a ab = diag(l, - 



° 2 (*) 
1 - Kr 2 



-a 2 (t)r 2 , -a 2 (t)r 2 sin 2 



(57) 

Note, that K = 0, ±1 is the spatial curvature for a flat, 
closed and open model respectively. This perfect dust- 
like configuration of the Riemannian limit can be trans- 
planted in a Finslerian set up by inserting the metric 
tensor (|57l) in relation (I55|) . Let us assume the com- 
mon Finsler case where the vector field y a is identified 



as the velocity of the congruence in question. Then, the 
4- velocity u a of the fluid flow stands for the anisotropic 
variables y a . By virtue of relation (|56[) the vector field 
b a is written into its irreducible parts 



Pua + b a , 



(58) 



where b a = h b b b stands for a preferred direction in the 
instantaneous rest space of the fundamental observer. 

In order to keep close to the FRW symmetry we impose 
homogeneity and isotropy. This, assumption switches off 
the spatial part of b a . Also, our analysis is restricted to 
monotonically decreasing functions of /3(a) to recover the 
FRW limit for large values of the scale factor. We con- 
sider the particular profile, j3 — cia~ n , where c\ stands 
for a small constant. However, the following analysis 
holds for a large class of profiles. Projecting along the 
4-velocity of the comoving observer relations (|55|) give 
back the evolution law for the scale factor 

la 2 -~K P a 2 (l + 4p)+Kp = -~K, (59) 

together with the linear independent continuity equation 



2/3 



1 + 2(5 



P = Q, 



(60) 



where /? remains an undetermined function of time. 

The Newtonian kinematical analogue of relation (|59|) 
is the motion of a point particle with kinetic energy 
K = id 2 , and total energy E = — \K. We can directly 
integrate the first order differential equation (|60|) to ob- 
tain the solution for the energy density 



P = Po 



1 + 2/3 



(61) 



In that case, our hypothetical particle falls in the poten- 
tial 



U=-\ Kpa a-^ + K P . 



(62) 



The modified FRW potential (|62l) corresponds to a col- 
lapsing dust-like fluid if the scale factor a is a decreasing 
function of time. On the other hand, it is natural to 
expect that the Finslerian contribution fades out as the 
fluid expands. Then, to recover the FRW behavior for 
large values of the scale factor we must impose that (3 
decrease faster than a -1 (n > 1 for our profile). In that 
case, for late times of the collapsing process the (3 con- 
tribution in the rhs of (1621) can play the role of a gravita- 
tional repeller. The main effect comes from the coupling 
term with the spatial curvature and depends on the Rie- 
mannian 3D-geometry (K = 0, 1, —1). In the same line 
with Raychaudhuri's equation (|45l) . we retrieve again a 
color-curvature contribution since /3 is the quantity that 
breaks the quadratic restriction (|56[) . Moreover, to re- 
sist gravitational pull of matter the spatial curvature K 
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FIG. 2: A sketch for the mechanical analogue of a test-particle 
moving in the potential (|62|) for a closed spatial geometry with 
P > and n > 1. The straight line represents the total energy, 
the dash line is the standard FRW potential and the solid one 
stands for the modified case (|62p . An initially contracting re- 
gion starting from the right of Umin will accelerate and then 
decelerate until the bouncing point. After the bounce, it en- 
ters an accelerating expanding phase until it crosses U m in- 
Then, the region decelerates until the self- gravity of the fluid 
dominates. The dust-ball recollapses and the process starts 
over revealing an oscillatory behavior for the 3D-volume ele- 
ment. 



and the Finslerian contribution /3 must be of the same 
sign. Notice, that the main effect fades out in case of flat 
spatial section. 

Consider a contracting region of our hypothetical 
space-time where initially the scale factor lies in the stan- 
dard FRW region of the potential (|62j) . Due to the self 
attraction of matter the contraction will be in an accel- 
erating phase as in standard gravitational physics. As 
the region shrinks, the coupling term between the spatial 
curvature and color starts to dominate and acts as an ef- 
fective pressure against the pull. Eventually, it forces the 
collapse to a decelerating phase until the fluid bounces 
to the potential at the turning point 



(ib 



np 



3(l + 2/3 & ) 



(63) 



where the subscript b denotes the value at the bounc- 
ing point. After the bounce the medium expands in an 
accelerating manner that turns to a decelerating phase 
until the gravity dominates again. Then, the fluid recol- 
lapse and the process repeats itself to infinity revealing 
an oscillatory behavior for the scale factor (see figure- 
[5]) • A similar behavior is retrieved in the context of very 
special relativity for the Bogoslovsky's line element [27j 
and on tangent Lorentz bundles [HJ]. The bounce is also 
possible in an open model if ft < but an expanding 
phase will never recollapse. Finally, in case of flat spatial 
sections the singular fate of a collapsing dust-like fluid is 
unavoidable. 



VIII. DISCUSSION AND CONCLUSIONS 

In summary, the intrinsic anisotropy imposed by the 
non-quadratic metric function F(x,y) gives birth to var- 
ious kinematical complications. The evolution of the de- 
viation vector along the medium's flow lines monitors 
the internal deformation and involves up to second or- 
der terms of £ a . Combined with the three distinct Ricci 
identities ([2"B)) - (|2"5j) that we get from the two types of 
covariant derivative, the evolution equations for the ex- 
pansion, shear and vorticity of the flow seems difficult to 
be given in closed form expressions. At this point one 
can either treat the tangent bundle as his base mani- 
fold and split his tensorial quantities to vertical and hor- 
izontal components (for textbook treatment see pj), or 
proceed to further geometric assumptions to get some 
analytical or numerical results. The first approach gives 
analogue relations with GR in the horizontal and vertical 
sub-bundles of TTM. providing a way to study focusing 
of congruences in the total space of the manifold (30j . 
However, one has to "translate" the covariant TTM- 
expressions back to the actual deformation of a medium's 
cross-section. 

On the other hand, at a first "brute-force" attempt we 
may try to solve the global problem for the congruence 
and its connecting vector for particular geometric con- 
ditions. Then, we can track back the deformation of a 
cross section through the evolution of the deviation vec- 
tor with respect to the arclength. The latter is given by 
the absolute differentiation along medium's flow lines and 
involves two distinct covariant derivatives. The first one 
is along the tangent to the flow lines u a = ^j— , while the 
other is along the evolution direction of the element of 
support, £ a — 21—. In the main bulk of the article we re- 
strict our analysis for simplicity to normalized time-like 
congruences that depend solely in the position coordi- 
nates, u a — u a (x). This particular class of Finslerian 
integrable curves is widely used to retrieve simpler for- 
mulas and isolate particular properties of the manifold 
[lj . They are closely related to the Riemannian ones and 
we can Lie drag a tensor field along them in a simple and 
straightforward way. However, the M a -bundle still "feels" 
the color of the structure mainly through the Cartan ten- 
sor. In fact, this is clearly depicted in the constrain re- 
lation for the Killing vectors (|3"3")l where the bundle can 
anisotropically expand in contrast to GR. 

The Lie derivative plays a central roll since we define 
the deviation vector as a Lie transported spatial congru- 
ence along the flow, £ u £a — 24]. Even in case of a 
deviation vector that depends solely in the position the 
evolution of a cross-section involves up to second order 
terms with respect to £ a . For a given direction l a , the 
first order deformation is given by the horizontal deriva- 
tive of u a and the second order by the Cartan tensor ([57)1 . 
The example of an isotropically evolving volume element 
is given where the expansion of the flow cannot be deter- 
mined ab initio. In a simplistic scenario we solve for the 
volume of the congruence and demonstrate how an expo- 
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nentially expanding region deflates or reach a caustic due 
to the Finslerian contribution. This suggests that emer- 
gent effects of QG may have a crucial kinematic reaction 
to inflationary models of the early universe and deserves 
further investigation. 

Possible GR modifications are more likely to appear in 
highly curved regions where QG physics may emerge. In 
GR, such conditions favor the formation of a singularity 
and/or a caustic of the congruence. Hence, the break- 
ing of the local symmetry of GR opens up the possibility 
to evade the focusing of neighboring particles and the 
space-time singularities (see for example [H, HH). From 
a kinematical perspective, any modification of the cur- 
vature theory will directly affect the deformation of the 
w a -bundle [22| . That is the case also in Finsler geometry 
where the curvature theory is extended with the intro- 
duction of color. The later monitors departures from the 
Lorentz symmetry since it reflects the breaking of the 
quadratic restriction on the distance module. Hence, in 
Finsler geometry the internal motion of a time-like flow 
is affected by the color (LV) and its variations. 

We concentrate our analysis in the expansion of a time- 
like flow. The complexity of the general problem forces us 
to pick a particular class of congruences along which the 
preferred direction imposed by the element of support is 
parallel displaced. Then, for a position dependent time- 
like bundle Raychaudhuri's formula can be retrieved. In 
that case, general conclusions for the expansion dynamics 
can be derived. The more interesting Finslerian effect is a 
coupling between curvature and color. This term is more 
likely to become dominant in high curvature environ- 
ments; exactly where we speculate that QG effects may 
become important. Under certain conditions it can repro- 
duce a repulsive effect to neighboring particles or add to 
the gravitational pull. The focusing theorem for our re- 
stricted congruence is given where irrotational geodesies 
of dust like matter can avoid the caustic. Nevertheless, 
the caustic prevention does not guarantee that the real 
space-time singularity vanish. However, it points out a 
possible QG mechanism that can slow down or even stop 
gravitational collapse. The formalism can be extended to 



study the effect on the shear and vorticity of the medium. 
In the same context the covariant perturbations over an 
almost FRW model can give back some imprints of color 
in the spatial distribution of the cosmological fluid 12311 . 
Also, the Raychaudhuri's equation for null geodesies [32[ 
would be desirable since it is of central importance in 
gravitational lensing and distance measurements [33j . 

One of the most famous solutions of Einstein's the- 
ory of gravity is the FRW space-time. It describes an 
isotropic and homogeneous fluid and it is a good approx- 
imation for the cosmological medium and for the interior 
of a star. In GR it clearly monitors the gravitational pull 
of matter through the evolution of the scale factor. The 
contracting phase is always accelerating and leads to the 
Big-Bang singularity, a — > 0. In phenomenological mod- 
els of quantum gravity and for non-standard matter fields 
the Big-Bang singularity is avoided in various examples 
34| . Usually, the medium before it reaches the singular 
point bounces back and enters an expanding phase. This 
models are candidate alternatives to inflation and they 
can also prevent the formation of the black hole horizon 
p35j . We recover the same behavior for an almost FRW 
model using the Randers metric in the osculating Rie- 
mannian limit. Moreover, the bounce occurs only when 
the space-time is spatially curved. The mechanism that 
prohibits the formation of the singularity is a coupling be- 
tween the spatial curvature and the variable that breaks 
the Lorentz symmetry. This result is in alliance with 
Raychaudhuri's equation and reveals from a geometric 
perspective that departures from the local symmetries of 
GR may lead to non-singular space-times. 
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